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^ ' Traffic congestion is usually observed at the upper stream of bottlenecks such as tunnels. 

Congestion appears as stop-and-go waves and high-density uniform flow. We perform simu- 
^SJ ' lations of traffic flow with a bottleneck using the coupled map optimal velocity model. The 

bottleneck is expressed as a road segment with speed reduction. The speed reduction in 
the bottleneck controls the emergence of stop-and-go waves. A phenomenological theory of 
bottleneck effects is constructed. 



KEYWORDS: traffic flow, bottlenecl<, optimal velocity, coupled map 



o 
o 

^^' 1. Introduction 

Traffic flow phenomena have been attracting scientific and engineering research interests 
\^ ' since the popularization of cars in the 1950s. The physical understanding of traffic flow on 

^N I expressways has been improved mainly on the basis of mathematical models and their com- 

<^ , puter simulations since the early 1990s. ^'^^ Many interesting features have been studied from 

f^ ■ the viewpoints of nonequilibrium statistical physics, pattern formation and transportation 

^ ' L 

^— V ■ phenomena. 

c/3 ' One interesting feature observed in traffic flow is the emergence of traffic congestion. Near 

o : 

•^ , a bottleneck, we observe high-density flow, which shows complex behavior in a density-flux 

^^- diagram. ^^ The high-density flow breaks down to stop-and-go waves at a distance from the 

M-c bottleneck. It is pointed out, however, that a bottleneck is not the origin of congestion. '^^ 

A bottleneck just increases the density of traffic flow. If the induced density is low, cars run 
smoothly. The uniform flow beyond a critical density is unstable and breaks down to stop- 
and-go waves. Sugiyama and Nakayama reproduced this feature using computer simulations. ^^ 
Mitarai and Nakanishi discussed that convective instability is closely related to the break-down 
to stop-and-go waves. ^"^' The increase in density due to the effect of a bottleneck is the key 
to understanding the emergence of congested flow at the upper stream of a bottleneck. In 
this study, we perform simulations of a system with a bottleneck under open boundaries and 
discuss the effect of the bottleneck on the emergence of congestion. 

Physical models of traffic flow are, in general, divided into two types; macroscopic and 
microscopic. The macroscopic models treat traffic flow as fluid. The microscopic models treat 
individual cars as particles and describe interactions among them. One of the microscopic 
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models is the Nagel-Schreckenberg model/'^' which is a cellular automaton model of traffic 
flow. Another type of microscopic model is the car-following model. ^-^ 

The optimal velocity (OV) model^^^ of traffic flow is one of the car-following models. The 
most important feature of the model is the introduction of the optimal velocity. In the OV 
model, each car controls its speed to fit the optimal velocity, which is decided by the headway 
distance to its preceding car. The model is described as a set of differential equations for the 
positions of cars. The model is suitable for treating the instability of the high-density traffic 
flow at the upper stream of a bottleneck such as a tunnel. ^^ 

We construct a simulation system for observing the emergence of congestion near a bottle- 
neck. The system should be an open-road system with the injection and ejection of cars. So, 
we employ the coupled map optimal velocity (CMOV) model of traffic flow,^^^ which is a tem- 
poral discretization of the OV model. The CMOV model is suitable for computer simulations 
with open boundaries. 

The organization of this paper is as follows: First, we describe the CMOV model and the 
setup of the simulations in §2. We construct a one-lane open-road system with a bottleneck. 
The bottleneck is implemented as a road segment with speed reduction and suppresses the 
flux in the bottleneck indirectly. The simulation results are shown in §3. Typically, stop-and- 
go waves are observed at a distant upper stream of the bottleneck. Near the bottleneck there 
is uniform traffic flow. We summarize the relation of the speed reduction in the bottleneck to 
the appearance of the stop-and-go waves. A phenomenological theory of bottleneck effects is 
discussed in §4. Section 5 is devoted to summary and discussion. 

2. Model and Simulation Setup 

We employ the CMOV traffic flow model, ^^^ which is a temporal discretization of the OV 
model. ^'^•* The CMOV model updates the position x{t) and the speed v{t) of a car by 

x{t + At) = x{t) + v{t)At, (1) 

v{t + At) =v{t) + a (^optimal (Ax) -v{t)) At, (2) 

where Ax is the headway distance to the preceding car, At is a discrete-time unit given as 
0.1 s in this paper, and q is a sensitivity constant. Each car controls its speed to fit the optimal 
velocity decided by the OV function 14ptimai(Ax), which depends on the headway distance 
Ax to the preceding car. The OV function is, in general, a sigmoidal function of the headway 
distance. For realistic simulations, we use the following form: 



Kjptimal (Ax) 



, , Ax — d , 

tanh ( 2 ) + c 

w 



(3) 



where parameters Vmax^d^w and c can be obtained through observations of the car- following 
behavior. We use the set of the parameters in Table I, which is compatible with that in ref. 
13. 
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parameter 


value 


unit 


d 


25.0 


m 


w 


23.3 


m 


'^max 


33.6 


m/s 


a 


2.0 


s-i 


C 


0.913 





Table I. Parameters in optimal velocity function eq. (3). 



L 



Injection l-^> 



J-^- Ejection 



Lb 
Bottleneck 

Fig. 1. Schematic view of system. Cars are injected from the left side and ejected away from the right 
side. A bottleneck region is located at the right end of the system. In our simulations, the length of 
the system and that of the bottleneck region are L = lOOOO(m) and Lb = 2000(m), respectively. 



Cars should stop to avoid backward motion and collision with preceding cars. The optimal 
velocity is negative if the headway Ax is less than Axmin which satisfies K)ptimai(^a;min) = 0. 
The avoidance is expressed by replacing eqs. (1) and (2) with 



x{t + At) = x{t), 
v{t + At) = 0, 



(4) 
(5) 



for Ax < AXinin- 

We construct a one-lane road of length L with open boundaries (Fig. 1). If a car arrives 
at the right end of the system, it is ejected from the system. The headway of a car following 
the car ejected from the system is set to be L as a headway sufficiently long. At the left end 
of the system, a car with zero velocity is injected if the distance between the left end of the 
system and the tail of the sequence of cars is larger than Axmin- 

We also introduce a bottleneck region of length Lb at the right end of the system, for 
observing its effect. The bottleneck is defined by reducing the maximum speed in the region. 
Namely, cars in the bottleneck run with the reduced OV function %^ijnar 



K,ptimal (A^) = rFopti^nal (Ax) , 

where r (0 < r < 1) is the degree of speed reduction in the bottleneck. 



(6) 
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3. Simulation Results 

We performed simulations with L = lOOOO(m) and L-q = 2000(m). After relaxation, we 
can see typical car trajectories in the space-time plane (Fig. 2). A high-density uniform region 
stably exists just before the bottleneck and maintains its length. The region is followed by 
striped patterns which correspond to stop-and-go waves. They propagate upstream, opposite 
to the direction of cars. No traffic jam emerges in the bottleneck. This feature was observed 
by simulation.^) Breakdown to stop-and-go waves was discussed as the convective instability 
of uniform flow without bottlenecks.^"^' 



Sr" 3200 9000 



1725 "^ ^ ..t^^ 
1800 , " 

\SQO) 







stop-and-go waves uniform flow bottleneck 

Fig. 2. Space-time plot of car trajectories with intermediate speed reduction r = 0.6. The horizontal 
axis denotes the positions of cars. The vertical axis denotes the time. The arrows represent regions 
of stop-and-go waves, uniform flow and the bottleneck, respectively. 



The inverse of the headway 1/Ax is plotted for each car as a snapshot in Fig. 3. The 
bottleneck induces three typical patterns of traffic flow. The emergence of stop-and-go waves 
depends on headways of cars just before entering the bottleneck. We can analytically discuss 
the linear stability of uniform flow using the OV model. ^^ The hatched areas in Fig. 3 show 
the density with which uniform flow is linearly unstable, 2^' jjjj^j^[(Aa;) > a. 

After relaxation for 2 h (72000 time steps), we calculate the average density at the 7800 m 
point by observing the flux and the average of the velocity per hour (36000 time steps) for 
each r value. The dependence of the average density pn on the speed reduction r is shown by 
the symbols D in Fig. 4. We define tl and r\] as the lower and upper bounds, respectively, 
between which pn remains in the hatched region. The boundary values tl and ry of the speed 
reduction are obtained as approximately 0.44 and 0.92, respectively, using the simulations. 
These values, however, are slightly different from the boundary values of the emergence of 
stop-and-go waves in the simulations. The discrepancy comes from the temporal discreteness 
of the CMOV model and the finiteness of the system length employed in this study. 

A test car is injected from the left end of the system to observe its behavior. Its typical 
trajectory in the head way- velocity plane is shown in Fig. 5 for intermediate speed reduction 
(rL < r = 0.6 < r\j). First, the trajectory draws a closed loop called the hysteresis loop while 
the car continues the stop-and-go motion during the approach to the bottleneck. As the car 
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Fig. 3. Snapshots of car density for r = 0.95 (a,), r — 0.6 (b) and r = 0.3 (c). The horizontal axes 
denote the positions of cars x (m). Each vertical axis denotes the inverse of headway l/Ax (1/kni). 
Each of the hatched areas corresponds to the headway with which uniform flow is linearly unstable. 
The uniform flow near the bottleneck indicated by an arrow, which corresponds to the uniform 
flow region in Fig. 2, breaks down to stop-and-go waves, as in case (b). 



approaches the uniform flow region before the bottleneck, the loop converges to a point on the 
curve of the OV function. Namely, the car moves with the optimal velocity given by the OV 
function. After the car enters the bottleneck, the trajectory moves to a point on the curve of 
the reduced OV function in the bottleneck. As the car approaches the right end of the system 
(the end of the bottleneck) , the uniformity of the headway in the bottleneck is lost (this part 
of the trajectory is not shown in Fig. 5). 

The relation between r and the density in the bottleneck, pB, is shown in Fig. 6. pB is 
observed at the 9000m point by the same method as pu- Except for the very weak speed 
reduction (r > r\j), pB is independent of r. We interpret these observed results in the density- 
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Fig. 4. Relation between speed reduction r and density pn (1/km) observed at 7800 m point. The 
observed values are shown as D. The hatched area corresponds to the density at which uniform 
flow is linearly unstable. The curve is given by our phenomenological theory discussed in §4. The 
speed reduction between tl — 0.44 and ru — 0.92 induces the density within the hatched region. 



flux relation. A uniform flow is observed in the bottleneck. We can calculate the flux q of the 
uniform flow with the optimal velocity in the bottleneck as a function of the density p: 



r ' optimal 



(7) 



Figure 7 shows the relation between q obtained from eq. (7) and the observed p^- pB corre- 
sponds to the maximum flux in the bottleneck. 

4. Phenomenological Theory of Bottleneck Effect 

We are interested in the effect of a bottleneck on the traffic flow at the upper stream 
of a bottleneck. Here, we construct a phenomenological theory of the bottleneck effect. We 
make two assumptions on the basis of our simulations. The first one is that a uniform-density 
flow exists just before the bottleneck. Thus, cars move with the optimal velocity just before 
entering the bottleneck. The flux giN entering the bottleneck is given by a function of pH'- 



Qmipu) 



PH Optimal ( 

\PH 



(8) 



where pu is the density of cars just before the bottleneck (Fig. 8). 

Second, after entering the bottleneck, pB is tuned to give the maximum flux: The flux in 
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Fig. 5. Typical motion of test car in plane of headway Ax(m) and velocity u(m/s) for case with 
intermediate speed reduction r = 0.6 (rL < r < ru). The car trajectory draws a hysteresis loop 
first. The trajectory converges to a point on the curve of the OV function T4ptimai(Ax) as the car 
approaches the bottleneck. 



the bottleneck qovT is given as a function of r by 



(b) 



gouT(r) = PbV^p^j^^^i 



1 

PB 



^^^PVopL.l[- 



PB = argmax pV^^l^^^ i^- 



(9) 



(10) 



The conservation law of flux requires that the flux qi^ entering the bottleneck is equal to 
QOVT- 

qmipn) = qo\jT{r). (11) 

Using this equation, we obtain pn just before the bottleneck as a function of r. 

The effect of the bottleneck is shown as the curve in Fig. 4 by solving eq. (11) numerically. 
The curve describes well the simulation results except for very weak speed reduction r > ru — 
0.92. From the curve, we obtain two boundary values r^ — 0.441 and r[j ~ 0.989 of the speed 
reduction. The intermediate speed reduction, r^ < r < r[j, induces the car density at which 
uniform flow is linearly unstable (the hatched area in Fig. 4). Thus, the stop-and-go waves 
emerge at a distant upper stream of the bottleneck. Therefore, by solving eq. (11) numerically, 
we can predict the occurrence of the stop-and-go waves using the value of r. 

The numerical value r^ of the lower bound agrees well with the simulation value tl- On 
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Fig. 6. Relation between density pB (l/km) in bottleneck (at QOOOm point) and speed reduction r. 
Except in very weak speed reduction cases r > r\j 2:^ 0.92, the density is independent of speed 
reduction. 



the other hand, the upper bound r[j disagrees with the simulation value r\j. The reason is 
explained simply. The injection method employed in this study cannot supply the maximum 
flux at the upper stream of the bottleneck. The maximum flux corresponding to the speed 
reduction for r > ry exceeds the injected flux. In other words, the flow injected from the left 
of the system is not sufficient to supply the maximum flux in the bottleneck. As a result, the 
assumption in the phenomenological theory is not satisfied for r > ru- 

5. Summary and Discussion 

We studied the effect of a bottleneck using simulations and a phenomenological theory. We 
employed the coupled map optimal velocity (CMOV) model for simulations. The bottleneck 
is defined as a road segment with speed reduction. We obtained the relation between the 
speed reduction r and the car density pu before the bottleneck (Fig. 4): The very weak speed 
reduction, r > ru — 0.92, does not increase the car density pu to form the stop-and-go waves 
(Fig. 3(a)). The very strong speed reduction, r < tl — 0.44, increases the density pu, which 
is high enough to stabilize the uniform flow (Fig. 3(c)). The bottleneck with the intermediate 
speed reduction, tl < r < ru, induces the high-density uniform flow (represented by the arrow 
in Fig. 3(b)) just before the bottleneck. This uniform flow is linearly unstable and breaks down 
to stop-and-go waves at the distant upper stream of the bottleneck. 

For the intermediate speed reduction, we find two important features. The first is that cars 
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Fig. 7. Curves denoting relation cq. (7) between flux q (l/5niin) and density p (1/kni) of uniform 
flow in bottleneck (at 9000ni point) for various values of speed reduction r. The curves correspond 
to r = 1.0, 0.9, 0.8, •• • ,0.2,0.1 (from top to bottom), respectively. The observed values of pB of 
the flow are shown as D. These values of pB correspond to those in Fig. 6. 
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Schematic diagram of effect of bottleneck. 



move with the optimal velocity just before the bottleneck. The second is that the bottleneck 
induces the maximal flux within the bottleneck itself. 

We employ these two features as assumptions for the phenomenological theory of the 
bottleneck effect. Using the conservation law of flux, we obtain the density just before the 
bottleneck as a function of speed reduction (the curve in Fig. 4). If the density corresponds to 
that of the linearly unstable uniform flow, the stop-and-go waves emerge at the distant upper 
stream of the bottleneck. Namely, we can predict the occurrence of the stop-and-go waves 
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using the speed reduction. 

The effects of bottlenecks have been studied previously.^' ^"^'^^•' The breakdown of high- 
density flow to stop-and-go waves has been observed in simulations.^' ^^' The breakdown effect 
has been discussed in relation to the convective instability of uniform flow.^"^-* The discussion 
in this paper, however, is based only on the relation between the speed reduction r and the 
bounds of the linear instability in ref. 4. The properties of the convective instability may 
affect the detailed properties induced by a bottleneck, including the stability of high-density 
uniform flow near the bottleneck. We will discuss these features elsewhere. 
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